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We investigate the transient evolution of finite-frequency current noise after the abrupt switching 
on of the tunneling coupling in two paradigmatic, exactly solvable models of mesoscopic physics: the 
resonant level model and the Majorana resonant level model, which emerges as an effective model 
for a Kondo quantum dot at the Toulouse point. We find a parameter window in which the transient 
noise can become negative, a property it shares with the transient current. However, in contrast to 
the transient current, which approaches the steady state exponentially fast, we observe an algebraic 
decay in time of the transient noise for a system at zero temperature. This behavior is dominant for 
characteristic parameter regimes in both models. At finite temperature the decay is altered from 
an algebraic to an exponential one with a damping constant proportional to temperature. 
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I. INTRODUCTION 



One central issue of mesoscopic physics focuses on 
the transport of charge carriers through nanometer- 
sized structures where quantum effects play an essen- 
tial role. In past decades, this research field expe- 
rienced a tremendous growth. Not only the electric 
current but also the shot noise, which is associated 
with the charge quantization of current carrying ex- 
citations, can reveal valuable information about their 
actual charge. For instance, the fractional charge of 
quasiparticlcs in the fractional quantum Hall regime 
or the charge of Cooper pairs in superconductors can 
be recovered in the Fano factor, which is the ratio 
of the shot noise to electric current. Nowadays, in 
addition to the measurement of current- voltage char- 
acteristics and noise (current auto-correlation func- 
tion), in many cases even higher cumulants in a non- 
equilibrium steady state situation can be accessed as 
wellii>2 The corresponding theoretical tool to gain the 
information about all current cumulants is referred to 
as full counting statistics (FCS) and was developed 
and successfully tested on many free as well as inter- 
acting systems during the last 20 years However, 
in preparative non-equilibrium, where certain param- 
eters are changed rapidly, only the current has been 
extensively addressed so far both experimentally and 
theoretically.—^ A notable exception has been work 
[H on transient current fluctuations at equal times. 
At the moment, much effort is invested to access the 
FCS in these situations, but has not been successful 
even for the simplest models available. Instead of fol- 
lowing this route, we directly calculate the transient 
finite-frequency noise for two exactly solvable models. 
In particular, we provide a comprehensive analysis of 
the zero-temperature case to extract the effects due 
to shot noise only. In addition, the influence of ther- 
mal fluctuations is addressed for one of the models. 
Moreover, our calculations may serve as a benchmark 
for various numerical simulational methods such as 
the density-matrix rcnormalization group (DMRG), 



the functional renormalization group (FRG), or the 
Monte Carlo technique, which have already been ap- 
plied to some models closely related to those to be 
treated below 

Our paper is structured as follows. In the next sec- 
tion we present our models and the observables of in- 
terest. Sections IIIII and IIVI are devoted to our main 
results, namely the analysis of transient noise in the 
respective models. Section [V] summarizes our find- 
ings. The Appendices offer details of several lengthy 
computations. 



II. MODELS AND OBSERVABLES 

The two models of interest are the resonant level 
model (RLM), which is equivalent to the non- 
interacting Anderson impurity model (AIM)^ and 
the Majorana resonant level model (MRLM), which 
corresponds to a special parameter constellation of 
the interacting resonant level model (IRLM) and can 
be mapped onto the Kondo model at the Toulouse 
point. 13 These two models are rare examples of exactly 
solvable systems in non-equilibrium ! 14 ' 15 We consider 
them in a two-terminal setup describing a quantum 
dot that consists of one single electronic level cou- 
pled to two electronic reservoirs at different chemi- 
cal potentials. In the RLM, the lead electrons are 
treated as non- interacting fermions (Fermi liquids), 
whereas in the MRLM, depending on the system re- 
alization, they are one-dimensional (ID) interacting 
fermions (Luttingcr liquids) and, in addition, perceive 
a Coulomb repulsion with an electron on the dot if one 
starts with a resonant level in a Luttingcr liquid,— 
and non-interacting Fermi liquids in the case of the 
Toulouse point of the Kondo model.— A typical re- 
alization of the RLM is a quantum dot on the basis 
of a semiconductor heterostructure in the regime in 
which electronic correlations on the dot are negligible. 
Alternatively, one can think of quantum dots in the 
deep Kondo limit, the transport properties of which 
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are dominated by the resonant level physics J^r— In 
contrast, the MRLM can be relevant in dot-lead struc- 
tures composed of single-wall carbon nanotubes, in 
which the conduction electrons are strongly correlated 
and form Luttinger liquids . 21 ' 22 

A. Resonant level model 

In this case the system under consideration is 
purely non- interacting, (i.e., charge carriers are non- 
interacting both in the leads and the dot region). 
Consequently, the spin degree of freedom is irrelevant 
and therefore we suppress the corresponding index of 
fermion operators. Moreover, we assume that the dot 
region can only be occupied by one single electron. 
In real setups, this is justified if the quantum dot is 
sufficiently small and the spin degeneracy of the en- 
ergy levels is lifted, for instance by applying a strong 
external magnetic field so that transport can occur 
effectively only through one level. The Hamiltonian 
then reads 

Hrlm = H + Hpj + Ht , (1) 

where Hq specifics the contribution of the free lead 
electrons 

Hq= 22 e k,acl a C k ,a , (2) 
k,a=L,R 

with Ck.a denoting the annihilation operator of an elec- 
tron in lead a with momentum k. The second ingre- 
dient is the dot Hamiltonian given by 

H D = AdU. (3) 

In addition, we consider tunneling processes between 
the leads and the dot region, represented by the cor- 
responding Hamiltonian 

H T = 7a(*)[Vi(a = 0)d + H.c] , (4) 

a=LM 

where d and i\) a are the annihilation operators of the 
dot and lead electrons, respectively. We define the 
operator of the total current through the constriction 

as 

f(t) _M^M, ( 5 ) 

where the operator for the current between an indi- 
vidual lead a and the dot is given by 

/ a (t) = i7«(*)[^(*)d(t)-H.c.]. (6) 

Anticipating the sudden switching of tunneling we 
consider later, we already included an explicitly time- 
dependent tunneling amplitude 7(f). For reasons 
of clarity, we always assume a symmetric coupling 
7l(£) = lR.(t) = l(t). The asymmetric case can be 
treated as well, but the main physical effects to be 
discussed in this paper are unaffected by its concrete 
choice. 



B. Majorana resonant level model 

We now turn to an extension of the RLM that ef- 
fectively describes interacting systems. One interest- 
ing realization of this model is an interacting reso- 
nant level sandwiched between two electrodes in the 
Luttinger liquid phase at the interaction parameter 
g = 1/2—^2- In addition it takes into account the 
Coulomb repulsion between the resonant level and the 
leads. The Hamiltonian is given by 

HmLM = Hk + Ht + He , (7) 

where 

H K = AdU + J2 H o& a ] (8) 

a=L,R 

is again the kinetic part describing the localized dot 
level and ID interacting fcrmions modeled by the Lut- 
tinger liquids, and 

Ht= J2 7a(t)ty£(* = 0)d+H.c.] (9) 

a=L,R 

is the usual tunneling part. The additional 
terml^M 

H c = \cdU ^ i>t(x = 0)ip a (x = 0) (10) 

a=L,R 

is responsible for the Coulomb repulsion. In a gen- 
eral non-equilibrium setting, this model has not been 
solved exactly so far. However, it has been shown that 
the special choice Ac = 2-7T (= 2nvp if the Fermi veloc- 
ity of the lead electrons vp 7^ 1) leads to a Hamiltonian 
quadratic in fermionic operators after some transfor- 
mation steps, namely bosonization followed by a uni- 
tary transformation and re-fermionizationJ^ The re- 
sulting model, which can be mapped onto the Kondo 
model at the Toulouse point 2 ^ is called the Majorana 
resonant level model (MRLM) and possesses an exact 
solution! 17 i 26 After the series of transformations men- 
tioned above its Hamiltonian can be written down in 
the following way 

-Hmrlm = H K [£, 77, a, b] + H' T [£, 77, a, b] , (11) 
where 

H K [£, 77, a, b] = iAab 
+ ijdx [rj(x)d x r](x) + £(x)d x £(x) + V£{x)r){x)] 

(12) 

governs the dynamics of the free lead Majorana fields 
T](x) and £(x) and local dot Majorana fermions a and 
b, which are related to the original dot operator by 
d = (a + ib)/ s/2, whereas 

H' T % 77, a, b] = -i [ 7+ &f (a = 0) - 7 -a^ = 0)] (13) 
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is an interaction term modeling couplings between 
lead and local dot Majorana fermions. Here, we in- 
troduced the coupling constants j± = ± Jr- We 
take our operator for the total current through the 
constriction in Majorana fermion representational 

i(t) = - i [ 7+ (t)b(t) v (t) + 7 _ (t)a(t)m , (14) 

with special emphasis on the time dependence of the 
tunneling coupling. For the rest of this paper, we 
also specialize to symmetric coupling in this case and 
therefore have j-(t) = and define 7+(f) = 7(t). 
It has to be noticed that the splitting of the current 
into left and right contributions is only reasonable in 
our derivation starting from the resonant tunneling 
setup between Luttinger liquids. In the Kondo pic- 
ture this is not meaningful since this model describes 
the scattering of conduction electrons off a local im- 
purity, in which tunneling between the electrodes is a 
single-stage process (electrode-electrode with a spin- 
flip of the impurity), while electron transmission in 
the Luttinger set-up is a two-stage process (electrode- 
dot-electrode). One further difference concerns the 
interpretation of the dot energy in the Kondo case as 
a local magnetic field. Thus the dot magnetization in 
the Kondo picture corresponds to the dot occupation 
in the Luttinger setup. In most of the calculations 
presented below only fields at the tunneling point are 
involved (which means x = 0), therefore we suppress 
the spatial coordinate. 



C. Noise and current fluctuations 

In contrast to the intuitive nature of current flow- 
ing through a conductor, one has a certain degree of 
freedom in defining the time-dependent current noise 
spectrum S(Q) in a full quantum treatment of a tran- 
sient problem. We use the following, rather general 
definition which is directly related to the conventional 
noise definition in the steady stated! 

s(n) = J d(t- t'^w-^stt, f) , (15) 

with the irreducible current-current correlation func- 
tion 

s(t,t') = (i(t)i(t'))-(i(t))(i(t')) , (i6) 

which quantifies the fluctuations accompanying the 
current flow. £ denotes the domain in the space of 
time differences t — t' in which information about the 
current correlations is available. In the most obvious 
case of stationary state S = (— oo, oo) and the current 
correlation function depends on t — t' only. Therefore 
the noise as defined in Eq. (|15[) is time independent. 
In general S(Q) is a time-dependent quantity though, 
as we shall see later. We can express Eq. (fll))) in terms 



of current cross correlators between different leads a 
and P 

s af} (t,t') = (i a (t)ip(t')) - (i a (t)) (itf)) (it) 

so that we obtain the decomposition 

S ( t > t ') = l E W)S afl W). (18) 

a,/3=L,R=± 

Throughout this article, we consider the sudden 
switching on of the tunneling of the form ^(t) = 7#(i), 
where 0(t) is the Heaviside step function. The substi- 
tution of new variables r = t — t' and T = t + t' 
effectively restricts the integration range from — T to 
T and finally leads to the transient noise formula (em- 
phasizing the explicit time dependence) 

5(0, T)= [ dre iUT 5(r,T), (19) 

J-T 

which has to be evaluated for our two cases. In a 
steady state, all Green's functions only exhibit a de- 
pendence on the time difference r. Thus, we can im- 
mediately carry out the r integration to access the 
stationary solution, which has to be equal to the tran- 
sient noise in the limit of infinite time T, 

S stat (n)= lim / + dre inT S(T,T). (20) 

T— ¥00 J _rp 

This relation serves as a consistency check of our re- 
sults. The unit of current noise is given by 7r 2 rGQ, 
where Go = 2e 2 /h is the conductance quantum and 
r is the hybridization constant expressible in terms 
of the tunneling amplitude 7 and the electronic den- 
sity of states of the leads po, which is assumed to be 
energy independent for the rest of this article. This 
seemingly crude approximation is often called the wide 
flat band limit. For the RLM, we take the conven- 
tion r = 27T/9o7 2 , whereas for the MRLM, we define 
r = 7 2 /2 using po = l/(27r) which is required by the 
transformation procedure. One particular advantage 
of the definition (fTD"]) is that it can be easily applied to 
the experimental data in the form of time-dependent 
current traces as presented in Ref. 0. Nonetheless, 
the solution of the transient problem as shown below 
can be very efficiently adopted to any other definition 
of the transient current as well. 



III. NOISE IN THE RLM 

A. Adiabatic noise and transient current 
evolution 

Before approaching the problem rigorously, we 
attempt an approximate calculation of the zero- 
temperature current noise by assuming that it follows 
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the transient current adiabatically. This ad hoc ap- 
proach can only work well when the corresponding 
switch-on time r sw is much larger than the typical time 
scale of the current evolution, which is proportional to 
1/r. Nonetheless, we would like to look into the sud- 

I 



from the transient current formula given in Ref. 
into the generalization of the Schottky formula^ for 
zero-temperature, zero- frequency (shot) noise in a 
steady state with an energy-dependent transmission 
coefficient ) 29 ' 30 which is nothing more than the second 
cumulant of the corresponding charge transfer prob- 
ability distribution. Then, we obtain the adiabatic 
noise evolution as 

/+V/2 j 
=-T(u,t) [l-7>,t)] . 
-V/2 2n 

(22) 

The function S a (t) is symmetric with respect to 
both voltage and dot level energy A and contains only 
the difference between the left and right Fermi func- 
tions. In general, the time dependence shows up oscil- 
latory behavior with frequencies V/2± A. On the con- 
trary, the envelope is exponential in time so that the 
stationary value is reached after a time proportional 
to 1/r. Of course, this might be just an artifact of our 
approximation. That is why we would like to attempt 
an exact analytical solution of the problem below. An 
important point is that for a large enough absolute 
value of the detuning |A|, the current noise according 
to our definition becomes negative, which is depicted 
in Fig. [TJ This peculiar feature, which persists for the 
transient case to be studied below, has not been re- 
ported in the literature so far. We briefly want to 
turn our attention to the total transient current which 
shares this property, illustrated in Fig. [5J It is due to 
the fact that the transmission coefficient 7~(w,£), in 
spite of being properly normalized, can become nega- 
tive. Although a net charge backflow at intermediate 
times seems to be counterintuitive at first sight, it can 
be made plausible since both Fermi levels appear to be 
almost at the same height in the case of strong detun- 
ing (i.e., when A represents the largest energy scale 
of all adjustable parameters). Of course, this prop- 
erty also applies individually to both the left and the 
right currents. Just after switching on of tunneling the 
electrons of both leads start to populate the initially 
empty dot and at the very beginning both In and 
II have the same sign. Due to the very high energy 
difference |A| on very short time scales an overpopu- 
lation occurs. After that the current signs change and 



den switching case r sw = to obtain a qualitative pic- 
ture of what might happen to the transient noise. To 
achieve our goal, we insert the effective time- dependent 
transmission coefficient for the initially empty dot 



(21) 
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FIG. 1. (Color online) Adiabatic noise evolution at fixed 
voltage V/Y = 1 and zero frequency Q/Y = for varying 
| A| /r = 0, 1, 2, 5, 10 (red solid, orange long-dashed, green 
short-dashed, blue dot-dashed, and black dotted curves). 



a negative net current can be observed for a rather 
short time interval. Negative transient current has al- 
ready been discussed by the authors of Refs. [jfl) and 
[32| , but in these works, it arises only if the bandwidth 
of the leads is small enough, whereas in our case, the 
bandwidth is taken to be infinite. Other adiabatic 
schemes have already been proposed previously (see, 
for example, Ref. [33| and references therein). 
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FIG. 2. (Color online) Transient current at voltage V/T = 
1 for various detunings | A| /Y = 0, 1, 2, 4, 8, 10 (red solid, 
orange long-dashed, green short-dashed, cyan dot-dashed, 
blue double-dot-dashed, and black dotted curves). 



7>,i) = 



r 2 - Te-" (r cos [(a; - A)^ - (a; - A) sin [(u - A)t] ) 

(w - A) 2 + r 2 
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B. Transient noise evolution 



We now want to study the transient behavior of 
current noise at finite frequency in its most general 
form. We compare our results with a steady state cal- 
culation at finite frequency (corroborated by an FCS 
calculation at zero frequency)^ In addition, we pro- 
vide compact formulas for various limiting cases at 



zero temperature. The method of choice is the non- 
equilibrium Keldysh Green's function technique as it 
provides an intuitive physical picture for every single 
constituent of relevant equations. As a cross check 
we performed the same computation using the func- 
tional integration technique and obtained precisely the 
same results. The substitution of our current operator 
Eq. (j6]) into Eq. (fTT|) and assigning times t and t' to 
different branches of the Keldysh contour, followed by 
the application of Wick's theorem yield o 35 ' 36 



J 



S a p(t + ,t'_) = 7(07(0 [G^ d + (t' 7 t)G% (t, t') + Gp + (t', t)G+ d ~ (t, t') 

'(t',t)G+-(t,t') - G- p + (t\t)G+-(t,t>)] (23) 



~ G da' 



r 



with the general definition of the Keldysh time- 
ordered Green's functions 



G (C (M') = -i(r c ^ ( w4(t')) 



and the definition of the S matrix 

S = T c c~'^c dtfl T(t) 



(24) 
(25) 

(26) 



treats both operators on equal footing. The average 
in Eq. ([24| is taken with respect to the coupled sys- 
tem, while the average in Eq. (|25p is performed with 
respect to the uncoupled one. The next, somewhat 
tedious task is to evaluate the various Green's func- 
tions. To achieve that, we make extensive use of the 
following general relation for the RLM case, obtained 
by expansion of the S matrix to first order and subse- 
quent re-exponentiation, 



where ( and C' specify the respective dot and lead 
operators. Here, we use a compact notation which 
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' OcCX 



,(t,f) 



dsj(s) 



9ll(t,s)GZ(s,t') 



(27) 



r 



where the upper indices indicate the branch of the 
Keldysh contour (— /+ for the forward/backward 
branch) and the lower ones specify the lead/dot op- 
erators. It proves to be advantageous to express all 
Green's functions in terms of the full dot Green's func- 
tion D(t,t') = Gdd(t,t') and the free lead Green's 
functions g aa ' (t, t')M- The Dyson equation for the full 
Keldysh dot Green's function in matrix form is 



where 03 = diag(l, —1) is the third Pauli matrix and, 
for later use, we defined the even/odd tunneling self- 
energy as 

E±(t,f ) = 7(*)7(*') [g L L(t,t')±g RR (t,t')] . (29) 
The free lead Green's function in Fouricr-Kcldysh 
space read a 14 ' 38 " — 



joa'(w) = 2TTip 5 aa > 



-1/2 n a 
-I 71c, -1/2 



(30) 



b{t,t') = A>(t,0 

+ I ' dt x I dt2Do(t,t 1 )a 3 t + (t 1 ,t 2 )a 3 D{t2,t'), (28) 



where n a (uj) = np(uj — fi a ) represents the Fermi-Dirac 
distribution function of the respective lead electrode a 
with chemical potential /Lt a , while the free dot Green's 
function is given by 
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D (t,t')=e 



-iA(t-t') 



-i[6(t - t')(l - n ) - 6(t' - t)n ] m 

-i(l - n ) — i[6>(t' - t)(l - n ) - 6(t - t')n 



(31) 



r 



where n denotes the initial population of the quan- 
tum dot. Using the above relations, we finally obtain 
the irreducible current-current correlation function 



S(t 1 t') = ^[S 1 (t,t') + S 2 (t 1 t')}, 
where we defined 

Si(M') = D-+(t',t)j:+-(t,t') 



(32) 



and 



S 2 (t,t') = -2 -Re 



+ Z+ + (t',t)D+-(t,t') (33) 



J dhD^t'^ZZ+itut) 



x J dt 2 D R {t,t 2 )T,Z + {t 2: t') 



(34) 



It has to be noted that this formula splits into two 
major parts. S±(t, t') involves the sums of Fermi func- 
tions and depends on the initial dot occupation while 
S 2 (t,t') contains the differences of Fermi functions 
and is insensitive to the initial preparation of the sys- 
tem. To apply this formula, we now have to compute 
the various dot Green's functions. Thereby, we use 
the relation D H = D h + D R — D A and the versions 
of the Langreth theorem^ for the greater and lesser 
Green's functions 



D+- = (1 + £>*£*)£>+- (1 + D A T, A ) 
D-+ = (1 + D R Y, R )Dn + (l + D A Y> A ) 



D R ^ 
D R E- 



-D A , 
(35) 



where integration over the internal time variables is 
implied. Depending on the initial dot occupation, one 
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of these expressions simplifies tremendously. For an 
initially empty dot D^ + {t,t') = 0, whereas for an 
initially occupied dot we have Dq (t, t') = 0. The 
retarded and advanced Green's functions were calcu- 
lated in earlier work a 41 ' 42 by solving the corresponding 
Dyson equation, thus we only provide the results for 
reference 

D R (t, t') = -i6(t - t')e- iA( *- t ' ) e- r (*- t ') (36) 

and D A {t,t') = [D R (t',t)]*. These functions are in- 
sensitive to the initial dot occupation, which is re- 
flected in the fact that they are solely dependent on 
time differences. In our noise calculations, we first 
evaluate the time integral to get a formula which ex- 
plicitly contains the Fermi functions and thus applies 
at arbitrary temperatures. We then restrict ourselves 
to zero temperature and give the quite lengthy result 
in Appendix [X] The energy integrals of the first part 
Si(0, T) then have — oo as a lower boundary owing to 
the wide flat band limit, whereas the corresponding in- 
tegrations in the second part S 2 (fl,T) are performed 
on compact supports. The complete finite tempera- 
ture result is provided in Appendix [Bj As an initial 
condition, we choose an empty dot. 



1. Steady state solution 

We want to check our results by calculating the 
steady state noise independently and comparing it 
later with the limit T — > oo of the transient noise. 
Taking advantage of time-translation invariance and 
transforming to Fourier space, we obtain the follow- 
ing analytical formula 



5^(0)=7 2 / ^[i3- + (a;)G^-(a; + f2) + G^(a;)£> + -(u, + 0) 



2tt 



G^MG+"( W +0) - G^ + ( W )G+-( U + fi)]. (37) 



r 



Unlike in the time-dependent case, the Green's func- 
tions of Eq. (|3"T|) are easily accessible and are obtained 
by inverting the corresponding Dyson equation in ma- 
trix form. Using another formalism, the steady state 
noise spectrum for the RLM has first been calculated 
by the authors of Refs. [Hf and J43|, which is in ex- 



cellent agreement with our result. We note in pass- 
ing that, in comparison to our graphs, the authors of 
the aforementioned references obtained mirrored noise 
spectra with respect to fi on account of their slightly 
different definition of the Fourier transformation. As 
an additional check, we then specialize to the case 
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= 0, which indeed yields the same stationary re- 
sult as an independent derivation from the cumulant 
generating function. 



2. Limiting cases 



continuous switching function /(0) = 0. For VL 
±oo, we have 



Urn S(n,T) = -, 
lim 5(O,T) = 0. 



(41) 
(42) 



For the zero-temperature shot noise, we give com- 
pact, analytical formulas for various limiting cases by 
holding all other quantities fixed. The only terms that 
contribute are those of Si(f2,T). For V — > ±oo, we 
obtain 



These are the usual limits of the unsymmetrized noise 
in the steady state. We note that the aforementioned 
cases are all independent of the initial dot occupation. 
On the contrary, for A — > ±oo, we have for an initially 
empty dot 



lim S(n,T) 

V—>-±oo 



r 
T 



(38) 



which is accompanied by the saturation of the total 
current through the constriction at high voltage, 



lim (/(t)>=±~. 

V— >±oo Z 



(39) 



These two limits do not display any time dependence. 
It should be mentioned that this is generally not ex- 
pected in a model with finite bandwidth e c , where the 
short time scale behavior of the transient current is 
dominated by oscillations with a period of l/e c Z For 
T — > 0, we have an offset 



lim S(Q,T) = - ( 

A — > — oo Z 

lim S(fl,T) = 0, 

A— >+oo 



-FT 



(43) 



whereas for an initially occupied dot, the limits are 
reversed. This remaining dynamics of noise is under- 
standable since, in the former limit, a tunneling pro- 
cess is allowed only for an initially empty dot that 
can be populated by one lead electron, whereas in the 
latter, an electron on the dot can jump to one of the 
leads. This jump probability is equal for electrons 
on/to both leads, thus the time-dependent net cur- 
rent vanishes although zero temperature fluctuations 
are present. All formulas are clearly in excellent agree- 
ment with the steady state result. 



lim s(n,T) 

T^0 



r 

4' 



(40) 



Long-time asymptotics: Zero temperature case 



This limit can be linked to the V — > ±oo case, which is 
the same, and could thus be interpreted as tunneling 
into vacuum. For an arbitrary switching procedure 
7(i) = J0(t)f(t), a detailed analysis shows that the 
offset is generated by a boundary term of the form 
oc /(0)/(T), which obviously disappears in case of a 



Apart from the special limits above, we now analyze 
the general long-time behavior of transient noise at 
zero temperature. The most astonishing feature is the 
temporal decay as a power law for large times. At zero 
frequency (f2 = 0), we obtain in case of an initially 
occupied dot 



' (47r) 2 ^, / „ (u- A) 2 + T 2 47r ^ J m (u - A) 2 + r 2 

v ; m,n=± J -°° y * m= ± J -°o 



T 2 f v/2 f v/2 (r 2 - - A)(a/ - A)) Tsinc [{u/ - u)T] ^ 



(2n) 2 J_ v/2 J_ v/2 [(c-A) 2 + r 2 ][(cy-A) 2 + r 2 ] 



(44) 



I 

where Si(x) is the sine-integral function, sinc(x) is the fies to produce 
cardinal sine function^ and g(T) comprises all terms 

which decay exponentially and are thus subleading in p2 i-o 1 + Si (wT) /(27r) 

T. For zero voltage at resonance (A = 0), this simpli- S(Cl = 0, T) = — / du; 5 — —5 h g(T). 

2tt J uj z + T- 

(45) 



To leading order in 1/T, we find that the transient 
noise evolution for large times is dominated by a power 
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law 



s(si = o,tt > i) 



(46) 



For increasing voltage, this distinctive feature grad- 
ually disappears until, at infinite voltage, we attain 
the limit of Eq. (J3HJ). It can only be retained by 
adjusting the detuning A in such a way that the 
Lorentzian peak in the integrand of the first term in 
Eq. (|4"4"| is shifted to the zero of one of the sine in- 
tegrals (i.e., to the position of one of the lead Fermi 
levels —V/2 or V/2). This tendency is illustrated in 
Fig. [3] Moreover, the feature is only dominant if the 
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FIG. 3. (Color online) Difference between transient noise 
and its steady state at V/T = A/T = Sl/T - 
for zero temperature (red solid curve) and finite inverse 
temperature /3F = 200, 100, 50, 20 (orange long-dashed, 
green short-dashed, cyan dot-dashed, and blue double-dot- 
dashed curves). We include the reference function 1/(tt 2 T) 
(black dotted curve). 

frequency fulfills the condition Sl/T <C 1, which can 
be seen in Fig. U where we depict the transient noise 
spectrum at different times. We note the pronounced 
discrepancy to the steady state noise spectrum around 
Sl/T = 0, an indicator of the algebraic decay. Apart 
from that region, the curves are almost indistinguish- 
able for 7T = 20 on the plotted scale. In Figs. 15171 
we display the effects of tuning the various parame- 
ters of the model, namely voltage, dot level energy 
and frequency for the case of an initially unoccupied 
dot. Obviously, one recognizes the gradual approach 
to the limits calculated before. We stress that, using 
our definition of noise, we still observe negative tran- 
sient noise in two important cases: large negative fre- 
quency or large positive/negative detuning for an ini- 
tially empty /occupied dot, although the steady state 
noise is always strictly positive. This is consistent with 
very small overall noise levels in the corresponding lim- 
iting cases (A — > ±oo and SI — > -co). Since at finite 
values of these parameters, shortly before approaching 
the extreme cases, we always have oscillatory behav- 
ior, we expect and indeed observe an undershooting 
below the zero line. We want to present evidence 
of a relation between the long-time asymptotics and 
a feature of the steady state solution. Indeed, it is 
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FIG. 4. (Color online) Transient noise at times Tr — 
2.5, 5, 10, 20 (red long-dashed, orange short-dashed, green 
dot-dashed, and blue dotted curves) and steady state noise 
(black solid curve) at V/T — 2A/F = 10 as a function of 
frequency Q/T. 
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FIG. 5. (Color online) Zero-frequency transient noise at 
fixed voltage V/T = 10 for various detunings A/T = 
—20,-10,-5,5,10,20 (red solid, orange long-dashed, 
green short-dashed, cyan dot-dashed, blue double-dot- 
dashed, and black dotted curves). Note the dominance 
of the algebraic decay of the green and cyan curves (V = 
±2 A). 
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FIG. 6. (Color online) Transient noise at resonance 
(A/r = 0) and fixed voltage V/T = 10 for various frequen- 
cies Sl/T — —20,-5,-2,2,5,20 (red solid, orange long- 
dashed, green short-dashed, cyan dot-dashed, blue double- 
dot-dashed, and black dotted curves). 



9 



0.20 

NO 

o 

tn 0.15 



0.05 



2 4 6 8 

T [r - '] 



FIG. 7. (Color online) Zero-frequency transient noise 
at resonance (A/F = 0) for various voltages |V|/r = 
1,2,5,10,20,40 (red solid, orange long-dashed, green 
short-dashed, cyan dot-dashed, blue double-dot-dashed, 
and black dotted curves). 



and the function r(T) incorporates all terms of sub- 
leading order (i.e., algebraic terms of higher order 
oc \/T a with a > 1 and exponentially decaying func- 
tions). Our conjecture Eq. ([47]) obviously reproduces 
our analytical result from Eq. (|46j) . Provided it is cor- 
rect, we conclude the dominance of the algebraic de- 
cay for such parameter constellations in which the dot 
level coincides with a Fermi level of the electrodes and 
its gradual disappearance for growing detuning of the 
dot level away from a Fermi edge. This is supported 
by our calculations as well as numerical evaluation, es- 
pecially by the limiting cases V — > ±oo and A — > ±oo, 
where the feature is absent. 

At this point, we would like to address the similar- 
ities and differences to the calculation from Ref. 0, 
which addresses transient equal time current-current 
fluctuations in an RLM setup. There the calculated 
quantity is 



striking that the algebraic decay of the transient noise 
is dominant at zero frequency, where the stationary 
noise spectrum is non-differentiable, its first deriva- 
tive having a discontinuity 6S' . Inspired by the plots 
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FIG. 8. (Color online) Difference between zero-frequency 
transient noise and its steady state value for parameter 
pairs (V/T, A/r) = (0.5, 0), (0, 0.5), (0.5, 0.5), (0.5, l) (red, 
orange, green, and blue dotted curves from top to bottom) 
compared with the respective reference curves calculated 
according to the function SS' /(iT) (black solid curves). 

of Fig. [H it is tempting to suggest the following gen- 
eralization of our transient noise formula to arbitrary 
values of V and A, 



S(n = 0, TV > 1) = S stat (f2 = 0) + — + r(T) 



ttT 



(47) 



where the discontinuity is given by 



SS' = lim 



dS" stat (ft) dS stat (n) 



n-s-o+ dfl n->o- <9f2 

r 2 \ 2-<5v,o ' (48) 



1 E 



27T ^± V( A + cr ^/ 2 ) 2 + r 



2 1 r2 



5(r = 0,T) 



dn s(n,T) 



(49) 



that is, Eq. ([16)) taken at t = t' . Moreover, the t of 
Ref. |§ is related to our parameter by t — T/2. The 
procedure presented there consists of taking a time- 
dependent bias voltage and assuming its dynamics to 
be sufficiently slow so that an adiabatic approxima- 
tion can be applied. On the contrary, in our case the 
tunneling coupling is switched on instantaneously and 
thus infinitely fast and anti-adiabatic. 



4-. Long-time asymptotics: Finite temperature case 

We now want to address the calculation of transient 
noise for finite temperature. The result obtained after 
a cumbersome calculation is provided in Appendix |5] 
We here concentrate on the salient feature which con- 
sists of a modification of the temporal decay compared 
to zero temperature, which is now exponential. In- 
deed, we observe that the presence of thermal fluc- 
tuations introduces a new energy scale to the prob- 
lem on which the new damping constant is linearly 
dependent. In Fig. [31 we contrast these two types 
of decay. We point out that in these plots, we have 
subtracted the respective steady state values due to 
thermal Jolmson-Nyquist noise. An estimation of the 
finite-temperature damping constant is provided by 
r' = tt//3 so that the envelope of the transient noise 
for large times is cut off by a function proportional 
to e -7rT /^, where j3 is the inverse temperature. For 
more details, see Appendix [Bj This behavior is not 
unexpected as the transition from algebraic decay at 
zero temperature to exponential decay at finite tem- 
perature is a quite general phenomenon, which occurs 
in various systems and is not restricted to temporal 
evolution. As an example, we cite the spatial decay 
of Fricdcl oscillations, which follows a similar pattern. 
Furthermore, our result has a dramatic consequence 
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for eventual numerical simulations, which depend sen- 
sitively on the approach to steady state. We thus con- 
clude that these should be performed at finite tem- 
perature to reduce computational effort. From an ex- 
perimental point of view, it should be an observable 
effect, at least at sufficiently low temperature where 
the Fermi functions are not much smeared out so that 
one can detect the decrease of the damping constant 
as a function of temperature in different parameter 
regimes. 



5. Correlation function for dot occupation 

We want to mention an interesting similarity to the 
Fourier transform of the correlation function for the 
dot occupation, 

F(Q,T)= / +T rf(i-i')e in(t -'' ) (n d (t)n d (t')>- (50) 

J-T 

In an analogous calculation as before it can be shown 
that this function already displays an algebraic long- 
time asymptotics. For the special case V = A = Q = 
0, we find to leading order 

F(ft,7T»l)«4f- ( 51 ) 



However, it has to be stated that the charge suscep- 
tibility x(Q, T) exhibits a purely exponential decay in 
time already at zero temperature as it is related to a 
retarded Green's function and thus involves a commu- 



tator. Its definition reads 

X(0,T)= [ +T d(t-t') X (t,t'), (52) 

J-T 

where x(t, i') is a retarded Green's function given by 

x (t,t') = i6(t-t')([n d (t),nd(t')}). (53) 

This behavior is not surprising though. The charge 
susceptibility represents the response function to ex- 
ternal fields. One particular realization of such fields is 
a finite voltage across the constriction. The response 
is then the current through the system which, as we 
know, has an exponential behavior. 



IV. NOISE IN THE MRLM 

We proceed along the lines of the previous Section 
to evaluate the transient behavior of current noise in 
the MRLM. 



A. Transient noise evolution 

The transient evolution of the current was calcu- 
lated in an earlier work.— We use the same formalism 
and define the Majorana Green's function according 
to the following prescription 

G cc (t,t') = -i(C(i)C'(*')> = -i{t(t)t'(t')s) o (54) 
with the usual definition of the S matrix 

£ = T c e^c d * 7(t)6( ^ (t) . (55) 

Hence, we obtain the irreducible current-current cor- 
relation function 



S(t + ,t'_) = \i(th(t') [G+-(t,t')G+ b -(t,t') 
1 



D+-{t,t')G+-{t,t') 



p(th(t>) [G^(t,t')G;+(t 7 t') D+-(t 7 t')g+-(t,t') 



(56) 



In the second line, we used the facts that the re- 
tarded mixed Green's function vanishes and that the 
77-Majoranas decouple from the transport process for 
symmetric coupling. For completeness, we write down 
the retarded Majorana dot Green's function which is 
obtained by solving the Dyson equation^! 



Df b (t,t') = -ie(t-t')f(t-t'), (57) 



where 

/(*) = ^H^' - r/ 2 )^' 4 + ( fi/ + r/2)e- n \ 

(58) 

with O' = ^(r/2) 2 - A 2 . We also use the Langreth 
formula Eq. ([35]) for D bb + {t,t') and express the mixed 
Green's functions in terms of the full dot Green's func- 
tions by an expansion of the S matrix and subse- 
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quent re-cxponcntiation. As a result, we find a similar 
structure of the irreducible current-current correlation 
function as in the RLM, 

S(t,t') = j[Si(t,t') + S 2 (t,t')} , (59) 

where we defined 

Si(M') = -D+-(t,t')E+-(t,t') (60) 

and 

s 2 (t,t') = J dt lJ D i? (t,t 1 )s:+(i 1 ,<') 

x J dt 2 EZ + (t,h)D A (t 2 ,t'). (61) 

The functions S± are defined in Fourier-Keldysh space 
as 



which is accompanied by the saturation of the total 
current through the constriction at high voltage, 



(67) 



lim (I(t)}=±-. 

V—>±oo 4 

At T — >• 0, we again have an offset 



lim S(n,T) = —. 



At this point, we would like to mention that the dis- 
crepancy of the latter result with Eq. (|4"0")) is due 
to a non-vanishing contribution from the first part 
of Eq. (|35p which is absent in the RLM case. For 
SI — > ±oo, we have 



and 



= l 2 g V r,(u) = 7 2 .9^(w) (62) 



E_(w) = 7 2 3«r,(w) = -7 2 5^(w) (63) 



where the free lead Majorana Green's functions are 
given b y 16 ' 17 



lim S(Q,T) = — (1 
n->+oo 4 

lim S(O.,T) = 0. 

Q— > — oo 



2e 



-rr\ 



(69) 
(70) 



However, for A — > ±oo, we have for an initially empty 
dot 



i / n' T . + n' - 1 



and 



(64) 



(65) 



The primes indicate that, instead of choosing the elec- 
trodes' real chemical potentials ^l,r = ±V/2, we have 
to insert effective ones fi' L R = ±V into the Fermi- 
Dirac distribution functions. As expected we recover 
the stationary state results of the authors of Refs. [13] 
and Gl. 



1. Limiting cases 

As for the RLM calculation, we give compact for- 
mulas for various limiting cases at zero temperature 
by holding all other quantities fixed. The contribu- 
tions are due to terms of T), again containing 
only sums of Fermi functions. In the following, we list 
them in the same order as before. For V — > ±oo, we 
obtain 



lim S(n,T) = -(l + 2e- 

V— >±oo o 



-FT 



(66) 



lim 5(fi,T) = ^e- 1 ?, 

A— y— co o 

lim S(n,T) = 0, 



(71) 



whereas for an initially occupied dot the limits are 
reversed. In relation to the RLM case, we state the 
qualitative difference that we have a temporal dynam- 
ics in the case of limits ft — > ±oo and V —> ±oo for the 
MRLM. We speculate that, at least in the IRLM case, 
the feature is due to the Coulomb interaction term in 
the Hamiltonian, which is absent in the RLM case. 
The seemingly slower exponential decay in the limits 
A — > ±oo with r/2 is not directly comparable to the 
RLM due to a different definition of T in both models. 
Of course, letting T — > oo in the above formulas, we 
find an approach to the expected steady state values. 



2. Long-time asymptotics: Zero temperature case 



In analogy to the RLM case, we identify a term with 
a similar structure involving sine integrals. For fi = 0, 
it is given by 
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S 1 (il = 0,T) 

r 2 



l " f +v , f +V , , (jfef2' + r/2)(mn' + r/2) r , A . 



64ft 

E 

cr=± 



<7=± 
V 

— oo 



8tt 2 



E 

m— ± 



; 2 (Si [(w - mV)T] - f 



r 



(w 2 + A 2 ) 2 + w 2 (r 2 -4A 2 



160'tt 



E 

m.ri— ± 



(0' - nr/2) 



im(ft + w) + (nQ' - T/2) 



ft(r), 

(72) 



where the function h(T) summarizes all terms that 
are exponentially decaying and thus subleading in T. 
Note that here, the voltage is doubled with respect 
to the RLM, a peculiarity due to the transformation 
steps from the original models. For V = A = 0, we 
come to the same conclusion apart from the prefactor 
and again obtain an algebraic decay, namely to leading 
order 



S{Q = o,rr> 1) 



1 



2tt 2 T' 



(73) 



In general, instead of one Lorentzian peak as for 
the RLM, the second term of Eq. (|72|) shows a two- 
peak structure with maxima at w = ±A. How- 
ever, this does not modify our conclusion. Obviously, 
the term has an appreciable effect only if V ps ±A 
(i.e., if the dot level almost coincides with one of 
the "dressed" lead Fermi levels). If the dot and Fermi 
levels move away from each other, the two peaks are 
no longer situated at the respective zeros of the sine 
integrals. For increasing f2, we also observe the grad- 
ual disappearance of this distinctive feature as shown 
in Fig. [9] Moreover, we emphasize that the transient 
noise as well as the current can also become negative 
in the MRLM. The transient noise evolution for vari- 
ous parameters in the case of an initially empty dot is 
shown in Figs. [T0lfT2l 
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FIG. 10. (Color online) Zero-frequency transient noise 
at resonance (A/T = 0) for various voltages \V\ /T = 
1,2,5,10 (red solid, green dashed, blue dot-dashed, and 
black dotted curves). 
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FIG. 11. (Color online) Zero- frequency transient noise 
at fixed voltage V/T = 10 for detunings A/T = 
— 10, —5, — 1, 1,5, 10 (red solid, orange long-dashed, green 
short-dashed, cyan dot-dashed, blue double-dot-dashed, 
and black dotted curves). Note the dominance of the al- 
gebraic decay of the red and black curves (V = ±A). The 
inset shows the same plot zoomed in the range [0,1]. 
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FIG. 9. (Color online) Transient noise at times TF = 
10, 20, 30 (red dashed, green dot-dashed, and blue dotted 
curves) and steady state noise (black solid curve) at V/T = 
A/T — 5 as a function of frequency Q/T. 



V. DISCUSSION AND OUTLOOK 

The most striking feature that distinguishes the 
zero temperature transient noise from the evolution 
of current and dot population in both the RLM and 
the MRLM is its algebraic temporal decay dominant 
for certain parameter sets. It achieves its maximum 
magnitude if one of the (dressed in the case of reso- 
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FIG. 12. (Color online) Transient noise at resonance 
(A/r = 0) and fixed voltage V/Y — 10 for various frequen- 
cies n/r = —20,-5,-2,2,5,20 (red solid, orange long- 
dashed, green short-dashed, cyan dot-dashed, blue double- 
dot-dashed, and black dotted curves). 



We expect this effect to survive in the case of re- 
alistic band structures beyond the adopted wide flat 
band limit since a finite bandwidth can only affect the 
transient behavior on short time scales. It is also in- 
dependent of the detailed switching mechanism as it 
is an effect at large times. However, we find that finite 
temperature destroys this effect by introducing a new 
energy scale determining the damping constant of the 
exponential decay. Thus, we expect the results to be 
observable at sufficiently low, but finite temperature. 

The possible avenues for further progress could be a 
detailed analysis of the impact of the on-dot interac- 
tions within the framework of the conventional Ander- 
son impurity model, or a discussion of transient noise 
in Kondo systems beyond the Toulouse point. 
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Appendix A: Analytical expression for transient 
noise in the RLM: Zero temperature case 

Below, we give an exact expression for the transient 
noise in the RLM at zero temperature for an initially 
unoccupied dot. The first part involves single integrals 
on the non-compact supports [-co, ±V/2] and reads 



with 



Sii ^ T)= ^ ± L dw U- A+ ^ 2+r2 5^ A ) 2 + r v ( ] 



si i(w,n,T) = — (r -e~ r7 Tcos[(A - lo - Vt)T] + e~ rT (A - u> - fi)sin[(A-a>- fl)T]) , 

47T 

s 1>a («, O, T) = ^ Yl 1 ■ ( Si ^ V / 2 + ") T 1 + Si ^ V / 2 + ") T 1 e " rT ) ' (A2) 

P>5=± 

ai , a (u,Sl,T) = V q-Si[((pV-q(Lu + A)-2n) + iqsT)T/2}c^- A ^ T / 2 . 

The second part consisting of double integrals, both on the compact support [-V/2, +V/2], is given by 

5 2 (0,T)= > / doj duj > — A3) 

irli-v/2 J -vii ~f K w - A) 2 + r 2 Uuj' - A) 2 + r 2 ] 



<T=± 
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with 

r 2 

s 2 i(u),u',n,T) = -— [r 2 - (u - A)(uj' - A)l Tsinc[(u/ - w- fl)T] , 
4tt^ 

r 2 

s 22 (w,u/,ft,T) = -e _rT [r 2 - (w- A)(w' - A)l cos [(w + u/ - 2A)T/2] Tsinc [(w' - w - 20)772] , 

r 2 

s 23 (w,a/,0,T) = +— e- TT (uj + uj' - 2A)Tsm[(uj + oj' - 2A)T/2] Tsinc \(u>' - u - 20)772], 
An 2 

s 2A (.M,n,T) - -11b. F- 1 "-?™:?-^ -T'- A) ^-«— ,-r„ _ 

47T 2 i(2w — uj' + A — 2S2er) + 1 



(A4) 

I 

Appendix B: Analytical expression for transient 
noise in the RLM: Finite temperature case 

In the case of finite temperature, everything is much 
more involved. We begin with the analog contribution 
to Si in case of zero temperature. It splits into two 
parts, one, where all integrals are performed, 



Si 



J] j 4i?r t Sl ( Z ' - <Bl ( 2 ' + 4i7re " 2rT [®1 (*, -XI*) - ®1 (*> 

"\ Too ^ \ „ — 2rTtvi /„ „ \1 /o:_ i ^ \ Too /_ _ „ — 2rToo /_ 



8tt 4/?r 

+ (2i7T - Xla) [*Bo (Z, —Xla) 



c~ 2rT <B 0, -ai,)] - (2i7r + xi*) [® n*) - e- 2rT <8 («, (Bl) 



4/3r 



1 ia; iCT 

7T — IW \ 

1 2 4tt 



IXlv 
47T 



and a second one where one frequency integration is 
left over, 

J 



Si 2 



r 2 



(T,T=± 

+ i*8 (z, -x 2<TT ) - iQS (2, x 2(TT ) + e 



n L M + Hr (w) I 
(2tt) 2 (w - A) + T 2 1 



-FT 



x 3o 



+ i*B (z, -a; 3crr ) - i*B (z, x 3aT ) 



e -rT/2 e -iT( w -A)/2 



27rn s f -7p) + iS (2, — £4o-r) - iS (z, a; 4crT ) 



(B2) 



e -rT/2 e iT(^-A)/2 



27ms 



X4 



iSo (-», -X4<Tr) - iS (z, X4o-t) 



Bi(«,o) B z {\ + with B de- 

noting the incomplete Beta function^ 
xi„ = p{(rV -2{A + Si)-2iT), x 2(TT 
£ (tV - 2fi - 2crw), x 3c7T = j3 (tV — 20 — 2aA), 



x 4aT — (3 (tV — 20 — aoj — a A — ioT), and 
z = c' 2 ^/ 13 . n s (x) = n F (-x//3) is the Sig- 
moid function^. The analog contribution to 5 2 is 
given by 
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r- 



Re y: 



duj til (u>) — (u) J /? 

(2tt) 2 i(w-A)+r 1 J/lar - J/i, T 



4— ^ ■ 

cr,r — ± 

e v «8 (z, y T ) + c v <8 (2, -y T ) + 2mn s {J^j - 2mn s (^J 



S (z, J/lo-r) - S (z, —yiar) 

(viaT ~Vt)T- 



2/3e" 



Z/2<tt - Vr 



Vr 



So (z, y r ) - e I *B (z, j/ 2<7T ) + V f - + J - ^ f - + g 



So (z, 2/ r ) + e 3 (-<B (z, -y T ) + 2i7m s 



(B3) 



where y T = P (tV - 2 A - 2iT), yi CTT = 
(rV - 2a; - 2erf}) , y 2CTT = P(rV-2u- 4crft) 
and 2/ 3(Tr = /9 (tV — w — A — 2crf7 - iT). To find the 
decay law of the noise correlation with time, one has 
to investigate term by term. First, we notice that all 
the remaining integrals are convergent even without 
the Beta functions (all the Beta functions are at most 
constant or decaying as a function of w). This is 
because of the overall Lorentzian-likc prcfactors. To 
estimate the asymptotics due to the Beta functions, 
the following power series representations turn out to 
be extremely useful (|z| < 1), 

B z (a, b) = z a V 9," b \ z" (B4) 
v ' ' ^ n (a + n) v ' 

n=0 v ; 



where (x) n = r (x + n) /T (x) is the Pochhammer 
symbol. In our case, we only need (1) = n! in the case 
of x = and (2) n = (n + 1)! in the case of x = —1. 
We introduce the following notation z' = f3 (r] u + i£T) 
where is a real function of oj and £ is a real constant. 
With g u we denote an arbitrary complex function of 
the variable u>. Then one obtains 



-tvT/B SIX. -"^ T 



-2-KnT/P 



g^e 7T1/P e 2 e 2 



5 n 1/2 + n - (Tp/ (Air) + i^p/ (4tt) 

1/2 + n - g£g/ (4tt) - (4tt) n _ 2mT/fi 



n=0 
oo 



n=0 
oo 



% (1/2 + n - (Tp/ (4^)) 2 + (77^/3/ (4^)) 2 



(B5) 



n=0 



Now, it is crucial that for every p there exists a posi- 
tive integer Np so that the modulus of the real part a n 
or the imaginary parts b n is majorized by 1 for n > Np. 
Hence, the real or imaginary part of the whole expres- 
sion can be estimated by a combination of finite poly- 
nomial p Np (e" 2 ^) and e - 2,r W / (l - c^ 2 ^). 

It is important to notice that the prefactor e^ rT//2 does 
not lead to an exponential increase in any case. The 
function g^ always suppresses this tendency. Although 
our result is valid for arbitrary temperature, we em- 



phasize that the zero temperature limit (i.e., P — > oo) 
is far from being trivial. 



Appendix C: Analytical expression for transient 
noise in the MRLM: Zero temperature case 

We provide the result for A ^ T/2 in the zero tem- 
perature case. Here, k = ± specifies the initially oc- 
cupied/empty dot. 
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2 poo paV < 

S x {£l,T) = Y J Si,i{Sl,T)+Y J dus lt3 (u;,n,T)+J2 dw£]a M (u,n,:r) (CI) 

uJaV i_J — oo 



i=l cr=± JcrV <r=±' 

with 

^ = 32^ £ _ 2A2 + 2 i.A^ (2 + Sl m + mV) + T] ) < 

m,ri— ± 



m.n— ± 



r 2 ^ 



si 5 (w, Q, T) = —j —5 - — -r- f Si [(w - mV - f2)T] - ^] 

8?r 2 (w 2 + A 2 ) 2 + w 2 (r 2 - 4A 2 ) V LV ; J 2/ 

si i6 (w,0,T) 



m— ± 

^2 



T 2 + nr/2)(m0 / + nT/2) ( ifc „-nn'-r/2)r/2 

32fi'2 7 r 2 + n (r/2 - ifcw)] [mft' + n (r/2 + ikw)] 



x (| + Si [(( 2f2 + p v - w ) + ifc + r/ 2 )) T /2]) , 



and 



r 2 r+V f+V 4 

^^^Mflv^ /_ v dw /_ v ^'E^c^^o,^ (C3) 



with 



fe ^ ± [° + k ( r / 2 - 1W )] i fi + ( r / 2 + 1W )] 



r 'f)Tl= V -™ (fT + nT/2) (m(]' + nT/2) (ifa- n tf-r/2)r/2 
2M ' ' ' j ^ [0' + n(r/2-iM] [mO' + n (r/2 + ifeu/)l 
x 2Tsinc [((2fi - 2ui' -uj)+\k (nO! + T/2)) T/2] , 
, ' o r r\ ST (n' + mT/2) 2 o ( i(u - u ')-2mn'-r)T/2 ( C4 ) 

x 2Tsinc [(2fi - w' - w)T/2] , 

S24 ( W) ^q,t)= V - -(n- + r/2)(^-r/2) (1(w ^)-r)T /9 

~± [0 + w (r/2 — iw)] [O' — m (r/2 + ia/)] 



X 



2Tsinc [(2fi - w - a/ + i2mfi') T/2] 
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